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Q i Abstract 

O ■ 

I ■ For superconducting coils with many turns, it is widely believed that 

5^ ' the magnetic field from magnetization currents is negligible compared to 

r^, the one generated by the whole coil. This article shows that this is not true 

^ • also for coils with a considerably large number of turns (up to thousands). 

We introduce a method to accurately calculate the AC loss in coils with 
C^ , virtually any number of turns that takes the magnetization currents into 

account. The current density and AC loss is analyzed for two windings 
made of coated conductor consisting of a stack of 32 pancake coils totaling 
' ^ . 768 and 6400 turns, respectively. 

lHj! 1 Introduction 

CN , The AC loss in ReBCCQl coated conductor coils has been a subject of intense 

> ■ study dlllllllElinilZllSlinilinilllllia- However most of the works are for single 

pancake coils or double pancakes, with the exception of [TOl[Tl]. This is because 
of the complexity of the calculations. Real windings may contain thousands of 

1/^ i turns, such as high- field magnets [13], SMES and transformers [2]. Therefore, 

it is needed to develop calculation methods for large number of turns without 

f~^ , degrading the accuracy. 

ff^ ' The present state of the art is to approximate that the effect of the whole 

coil in a certain turn is the same as an applied magnetic field. This applied 
field ("background magnetic field") is computed by assuming that the current 
density is uniform in the rest of the turns (we name this approach as "uniform 
approximation"). Afterwards, the AC loss in the turn of study is estimated 
j^ , by either measurements in a single tape [15j [16] or by numerical calculations 

[TT] . The problem of this approximation is that the neighboring turns shield the 
background magnetic field, in a similar way as in a stack of tapes [IHl [H]. It is 
a question how important is this effect in coated-conductor windings consisting 
on stacks of pancake coils, where the shielding effect is expected to be high. 

An important step forward for single pancake coils has been the continuous 
approximation [9j . This approximation substitutes the detailed coil cross-section 



^_ReBa2Cu307_a; where Re is a rare earth, usually Y, Gd, Sm or a combination of those. 



by a continuous one with a Jc corresponding to the cross-section average of the 
original coil. The current constrain is set by fixing the line integral in the 
axial direction to a certain value for every radial position. The continuous 
approximation is neither applicable to coils made of tapes or wires with non- 
negligible thickness of the superconducting core (such as Bi-2223 tapes or Bi- 
2212 or MgB2 wires) nor for pancake coils with a large separation between 
turns. 

In this article, we introduce a method to simplify the AC loss calculation 
without practically degrading the accuracy for general coils, including stacks of 
pancake coils with many pancakes, solenoids, and coils made of conductors with 
arbitrary thickness. In addition, we present the current density and the AC loss 
characteristics in two windings made of 32 pancake coils, with a total of 768 
and 6400 turns, respectively. 

The structure of the article is the following. First, we outline the general 
numerical method that takes into account the interaction of the magnetization 
currents between all the turns, which serves to monitor the accuracy of the 
considered approximations. Later, we introduce the simplified method. Subse- 
quently, we discuss the calculated current distribution and AC loss for the two 
example coils made of coated conductor and evaluate the error committed by 
the considered approximations. Finally, we present our conclusions. 

2 Calculation method 

In this section, we first outline the general numerical method, which calculates 
the actual current distribution in all the turns of the coil, and afterwards we 
present a method to speed up the calculations for large number of turns. 

2.1 General numerical method 

In this work, the basic numerical method to calculate the current density and 
the AC loss in the superconductor is the Minimum Magnetic Energy Variation 
(MMEV) method, as described in [5ni[SJ[S]. The results of this method agree 
with experiments for coils [5J [TU] . Additionally, MMEV has also been applied 
to model levitation [3T1 [52] . The MMEV method calculates the detailed current 
density, J, in the superconductor and assumes the critical-state model in its 
general form. That is, \J\ < Jc for E = and \J\ = Jc for \E\ > 0, where Jc is 
the critical current density and E is the electric field. In this work, we take the 
cylindrical symmetry of the coils. Although this method is much faster than 
commercial Finite Element Methods (FEM) [23] j calculations for coils with a 
large number of turns are still lengthy. 

2.2 Approximations for large number of turns 

Next, we present a model to speed up the calculations practically without de- 
grading the accuracy that is applicable to any numerical method, not only for 
MMEV but also for models based on FEM. 

In order to calculate the current distriution in a certain turn of the coil, the 
most straighforward approximation is to assume that the current density in all 
the other turns is uniform ("uniform" approximation). This approach neglects 
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Figure 1: Normalized AC loss {2ttQ/ jiqI^, where Q is the AC loss per cycle 
and tape length and Ic is the critical current) for one of the 200-turn pancake 
coils composing the "continuous" coil (table [Ij . The inset shows the distortion 
of the magnetic field generated by a single tape with constant Jc transporting 
20 % of the critical current relative to the magnetic field assuming uniform 
current density. The two curves are for lines parallel and perpendicular to the 
tape surface, respectively, with origin at the center of the tape (only the region 
outside the superconductor is shown). 



the magnetic shielding effect originated by the magnetization currents in the 
neighbouring turns, which often introduces significant errors. However, this is 
the approximation that is most commonly used in practice. 

A more accurate approach is to assume uniform current distribution only 
for the turns farther than the turn of study by a certain distance. We name 
this approach "neighbor approximation" . This approximation is based on the 
fact that the distance from the tape where the magnetic field is influenced by 
the non-uniformity of J is of the order of the width of the tape (insert of figure 
[T|) . For stacks of pancake coilf|j, it may be convenient to take each pancake coil 
as an indivisible unity. That is, J is either assumed uniform in all the turns of 
the pancake coil or the numerical routine calculates the detailed J in the whole 
pancake. The convenience of this approach is clear for coils that their radial 
thickness is smaller than the tape width. For closely packed pancake coils with a 
thickness larger than the tape width, it is useful to keep the pancake coils as an 
indivisible unity and assume the continuous approximation [5] for each pancake 
coil. Summarizing, for stacks of pancake coils, the neighbor approximation 
consists on assuming uniform J in all the turns except in the pancake where J 
is calculated and its neighbor pancakes up to a certain order (for instance 1st 
or 2nd neighbors). 

3 Studied coils 

In this article, we analyze two examples: one taking the real geometry of the 
pancake coils (real thickness of the superconducting layer and spacing between 
turns) and the other one taking the continuous approximation for each pan- 
cake coil. We name these coils "detailed" and "continuous" , respectively. Their 
number of turns and assumptions for J^ are outlined in table [T] For both coils, 
the width and thickness of the tape are 3.96 mm and 1.4 /xm, respectiveljo, the 
separation between turns is 188 and 465 /xm in the radial and vertical direc- 
tions, respectively, and the inner radius is 29.5 mm. For the "detailed" coil we 
take the anisotropy and magnetic field dependence of Jc from [10], related to a 
particular SuperPower tape [24] at 77 K. For simplicity, we assume a constant 
Jc for the "continuous" coil corresponding to 100 A of critical current (for the 
'continuous' coil, Ic = 100 A requires temperatures below 77 K). With this pa- 
rameters, the maximum magnetic field at the center of the coils is 0.3 and 4.7 T, 
respectively. The computations use 1 current element in the tape thickness for 
the "detailed" geometry and 20 elements in the radial direction of the pancakes 
for the "continuous" case, between 76 and 154 elements across the tape width 
(the highest value is for the lowest current amplitudes, /„), and a tolerance of 
Jc between 0.20 and 0.035 % of Jc{B = 0) (the lowest values are for the lowest 

4 Results and Discussion 

This section presents the results of the current density and AC loss for the 
two studied coils: the "detailed" and "continuous" coils (see table [Tj and sec- 

^ Solenoids can be approximated as a stack of pancake coils. 

^With the term 'tape' we mean the superconducting layer of the coated conductor, not the 
whole coated conductor. 



Table 1: The studied coils. The continuous approximation j9j is taken for the 
"continuous" coil. More data in the text. 

Coil name Turns Pancakes Assumption 

per pancake for Jc 

"detailed" 24 32 Dependent on B 

and its orientation, 
"continuous" 200 32 Constant 



tion [3] for more details) . The applicability of the "uniform" and "neighbour" 
approximations are also discussed. 

4.1 "Detailed" coil 

The current density in the "detailed" coil presents the following features (fig- 
ure [21 shows the average current density across the superconductor thickness 
for the central turn of each pancake: 12*'' turn from the inner radius). The 
uniform approximation fails to predict the current distribution (figure [2K)- The 
reason is that in the real case, the neighboring turns in the radial direction 
shield the magnetic field created by the rest of the turns, while the uniform 
approximation neglects this effect. As a result, for the uniform approximation 
the magnetization currents are larger and the plateau at the turn center disap- 
pears. In addition, the actual current distribution in the neighboring pancakes 
influences the magnetic field. In the real case, the neighboring pancakes par- 
tially expel the magnetic field, which concentrates between the pancakes. This 
increases the magnetic field close to the tape edges, locally reducing Jc- As a 
consequence, Jc decreases compared to Jc for the uniform approximation (see 
figure [2^) . On the contrary, taking into account the neighbor approximation of 
1st order already produces a satisfactory result (figure ^Bp) . The approximation 
improves with approaching to the coil ends. The cause is that the average dis- 
tance from the rest of the superconducting turns is the largest, and thence the 
magnetic field distortion from the J non-uniformity is the lowest. 

The validity of the uniform and neighbor approximations can be quantified 
by the error in the total AC loss (figure [S]). The uniform approximation largely 
over-estimates the AC loss because it over-considers the magnetization currents. 
The difference increases with decreasing the current amplitude. The coil rated 
current (critical current of the weakest turn) is computed to be 51 A, compared 
to 128 A for an isolated tape. Again, the neighbor approximation provides a 
satisfactory result. The obtained values are better than for the uniform ap- 
proximation, also for the neighbor approximation of 0th order. That is, J is 
assumed uniform in the tapes in all pancakes, except the one where J is calcu- 
lated. The neighbor approximation slightly under-estimates the AC loss. The 
error increases with decreasing the current amplitude because at low amplitudes 
the current concentrates close to the tape edges, and thence it is less uniform. 

4.2 "Continuous" coil 

The results for the "continuous" coil qualitatively describe real-size magnets. 
These results are the following. 
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Figure 2: Average current density across the tape thickness for the 12 turn 
starting from the inner radius for the "detailed" coil (table [IJ . The uniform 
approximation (unif. appr.) fails to predict the current density (a), while the 
neighbor approximation (neigh, appr.) is satisfactory (b). The calculated case 
is for the peak of the AC current with 46 A amplitude. Only the lower half of 
the coil is shown, although the whole coil is calculated. 
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Figure 3: The AC loss for the "detailed" coil shows that the results for the 
"neighbor" approximation converge to those for the "real" geometry with in- 
creasing the order of the neighbors involved, while the "uniform" approximation 
introduces a large error. 



The current density by assuming only the continuous approximation is in 
figure m Note that the continuous approximation practically does not introduce 
any error (see the AC loss for one isolated pancake in figure [T]) . Back to figure 
m the four top pancakes are saturated with magnetization currents, while in the 
other pancakes there appears a sub-critical region with roughly uniform current 
density. For the unsaturated pancakes, the critical region is wider at the mid 
radius, its width decreases with approaching to the inner and outer radius and 
sharply increases at those radius. The causes are that, first, the radial magnetic 
field presents a maximum close to the mid radius and, second, the magnetic 
shielding effect from the magnetization currents drops at both edge radius (see 
figure 4 in [TS]). 

The AC loss in the "continuous" coil without making further approximations 
presents a change in slope, from around 2.6 at low current amplitudes to around 
1.7 at high amplitudes (see figure [SJ. The current amplitude where this change 
of slope appears corresponds to the saturation by magnetization currents in the 
end pancakes, which contribute the most to the AC loss. Therefore, this change 
of slope evidencies that the magnetization loss is dominant. 

One may expect that the uniform approximation is accurate for magnets 
with many turns because in the turns that contribute the most to the AC loss 
(top and bottom pancakes (TU]), the local magnetic field is much larger than 
the self-field of one turn (for the situation in figure SI in the central turn of 
the top pancake the section average self-field is around 10 mT, while the radial 
component of the background magnetic field is around 530 mT). However, the 
uniform approximation largely over-estimates the AC loss except close to the 
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Figure 4: Upper half of the cross-section of the "continuous" coil (the model 
calculates the whole coil). The current density is for the peak of the AC cycle 
with amplitude 0.32/c- 
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Figure 5: The same as figure [3] but for tfie "continuous" coil (the "real" case 
takes only the continuous approximation[9]). The AC loss normalization is the 
same as in figure [1] 



critical current (see figure [S]) . The reason is that the uniform approximation 
neglects the shielding currents from the whole pancake, which are much larger 
than for one single tape. These shielding currents decrease the local radial 
magnetic field in the top pancake from around 530 mT for uniform currents to 
280 mT in the real case (at the boundary between J — +Jc and J — —Jc, figure 
S]). The effect of the shielding currents decreases with increasing the transport 
current amplitude. The cause is that with increasing the net current both the 
magnetic field created by the other pancakes increases and the shielding currents 
decrease. 

The error caused by taking the neighbor approximation is low, decreasing 
with increasing the order of the neighbors (figure [5]). However, this approxima- 
tion is less effective for the "continuous" coil as for the "detailed" one (compare 
figures [3] and [5]) because the pancakes for the former are one order of magnitude 
wider than for the latter, and thence the magnetic field due to the non-uniform 
current density expands to larger distances. 



5 Summary and conclusions 

Summarizing, in this article we have discussed the current density and AC loss 
in coils with a high number of turns, in particular two coils with 768 and 6400 
turns, respectively. We have also presented the neighbor approximation and we 
have found that the error that this approximation introduces is small. However, 
the usual uniform approximation (J assumed uniform except in the turn of 
study) fails to predict the current density and the AC loss. The reason is that 
the magnetic field generated by the magnetization currents is comparable to the 
background magnetic field, also for coils with thousands of turns. 



In conclusion, we have found that for coils made of stacks of pancakes with 
closely packed thin conductors, it is necessary to use relatively detailed cal- 
culations in order to achieve a satisfactory AC loss estimate. Combining the 
neighbor approximation with the continuous approximation, it is feasible to 
accurately calculate the AC loss in coils with virtually any number of turns. 
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